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' In a recent work, we presented the first application of tfie Poisson- Wiseman- Anderson 

metfiod of 'matched expansions' to compute the self-force acting on a point particle mov- 
ing in a curved spacetime. The method employs two expansions for the Green function 
which are respectively valid in the 'quasilocal' and 'distant past' regimes, and which may 
be matched together within the normal neighbourhood. In this article, we introduce the 
method of matched expansions and discuss transport equation methods for the calcu- 
lation of the Green function in the quasilocal region. These methods allow the Green 
function to be evaluated throughout the normal neighborhood and are also relevant to 
X * I ' a broad range of problems from radiation reaction to quantum field theory in curved 

{^j^^ spacetime and quantum gravity. 
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, 1. Matched Expansion Approach to the Self- Force 

'sj" i A charged particle moving in a curved background spacetime will not exactly follow 

C5 I a geodesic of that spacetime. Within the framework of perturbation theory, the par- 

' tide is seen at linear order to exert a force on itself, a self-force, which is responsible 

for the particle moving away from geodesic motion. 

A formal expression for this self-force has been derived^^* in terms of some easily 
calculated local terms and a tail integral over the entire past-history of the particle. 
For the case of a scalar charge, this tail term is given by 



o 
o 



X 

C^' 9'/ V^G,et(2(r),z(T'))dT', (1) 

where the world-line of the particle is given by z(t), q is the charge and G-[ct{x, x') 
is the retarded Green function. For electromagnetic and gravitational charges, the 
tail term involves a similar integral of the Green function along the past worldline. 

The greatest difhculty in computing ([l} arises from the computation of the 
Green function for points arbitrarily separated. For x' x, the Green function 
has a troublesome J-divergence behavior, while for x and x' separated by large 
distances, certain analytic expressions for the Green function are not valid. Within 
the matched expansion approach,^ these difficulties are avoided by splitting the 
integral ([T]) in two parts at a matching point and tackling each part independently. 
We have given an example calculation of the distant past contribution for the Nariai 
spacetime.^ We focus here on the contribution from the recent past, the quasilocal 
contribution. 
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Provided the quasilocal region is sufficiently small (within a normal neighbor- 
hood), the retarded Green function is given by the Hadamard form, 



G,ctix,x') =e^{x,x'){U{x,x')S{a{x,x')) -V{x,x')e{-<7{x,x'))} , (2) 



where U{x,x') and V{x,x') are regular, geometrical bitensors, 6 is the Dirac-delta, 
0_ is the step function and ensures causality and <t{x, x') is the Synge world- 
function. Noting that a ^ for all points in the normal neighborhood contribution 
to ([1]), it is clear that only V{x, x') will contribute to the quasilocal part of ([T]). 

2. Coordinate Expansion of the Green Function 

The bitensor V{x,x') may be written as an expansion in the coordinate separation 
of X and x' . For static, spherically symmetric spacetimes, this takes the form 



with 7 being the angular separation of the points. The coefficients Vijk{x) are calcu- 
lated using methods which rely on the symmetry of the spacetime.^ We developed a 
Mathematica code capable of calculating these coefficients to very high order (up to 
(Aa;")^°) and applied the theory of Pade approximants to improve convergence and 
extend the domain of validity of the series.* This allowed us to accurately calculate 
(Fig. [J) the Green throughout and up to the edge of the normal neighborhood, 
beyond which the Hadamard form can no longer be assumed to hold. 



Fig. 1. (Colour online) Left: The retarded Green function for points separated along the time 
direction in the Nariai spacetime. The regular Taylor series (blue dashed line) diverges before 
the edge of the normal neighborhood, which is at t — t' = 6.5. The Pade resummed series (solid 
red line) remains valid and closely matches the exact value (black dots) up to the edge of the 
normal neighborhood. Right: V{x, x') for a massless, scalar field along the light-cone of a point 
at r = lOM in Schwarzs child Spacetime. 



3. Covariant Methods for Calculating the Green Function 

Covariant methods for calculating the Green function of the wave operator and the 
corresponding heat kernel, are central to a broad range of problems from radia- 
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tion reaction to quantum field theory in curved spacetime and quantum gravity. 
Traditionally, most attention has focused on the diagonal value of the heat kernel 
K^Aix,x;s), corresponding to the coincidence limit x' ^ x of the Green function. 
By contrast, for the quasilocal contribution to the matched expansion self-force, we 
seek expansions valid for x and x' as far apart as geometrical methods permit. 

An extremely elegant, non-recursive approach to the calculation of the Green 
function has been given by Avramidi,^ but as his motivation was to study the 
effective action in quantum gravity he was primarily interested in the coincidence 
limit. We presented^*' Avramidi's approach in the language of transport equations 
and showed that it is ideal for numerical and symbolic computation. In so doing we 
built on the work of Decanini and Folacci^^ who wrote many of the equations we 
presented and implemented them explicitly by hand. However, calculations by hand 
are quite impractical for the very high order expansions we would like for self-force 
calculations. Instead, we used the transport equations as the basis for Mathematica 
code for algebraic calculations^^ and C code for numerical calculations.'^'^ 

Using our numerical code, we calculated V{x,x') on the past light-cone of a 
point X in the Schwarzschild spacetime. With V{x, x') a solution of the homogeneous 
wave equation, nV{x,x') — 0, this light-cone data could serve as initial data for 
a characteristic evolution of V{x,x') within the light-cone. For points within the 
normal neighborhood of x, there would be no caustics ~ a frequent cause of difficulty 
in characteristic evolution codes. Nevertheless, the allure of evolving past caustics 
and beyond the normal neighborhood is appealing. It may reveal new insight into 
the relation between the Hadamard form and the structure of the Green function 
outside the normal neighborhood. 
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